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Abstract 

We briefly discuss the status of the NNLL QCD calculations in the inclusive rare 
B decay B —> X s i + i~ . Two important ingredients, the two-loop matrix elements 
of the four quark operator (P 2 and the bremsstrahlung contributions, were quite 
recently finalised. The new contributions significantly improve the sensitivity of 
the inclusive decay B —* X s l + l~ decay in testing extensions of the standard model 
in the sector of flavour dynamics; for instance the two-loop calculation cuts the 
low-scale uncertainty in half and the bremsstrahlung calculation leads to a 10% 
shift of the position of the zero of the forward-backward asymmetry, 
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Figure 1: The decay -B —> X s l + l~ at one-loop. 

1 The Decay B — Ayr/r 

Inclusive rare H decays like I? —*• or B —> X s £ + £~ are very important tools to understand 

the nature of physics beyond the Standard Model (SM). The stringent bounds obtained from 
B —► Ay/ on various non-standard scenarios (see e.g. 10 0 ) are a clear example of the 
importance of theoretically clean FCNC observables in discriminating new-physics models. 

In comparison with the B —> A s 7 , the inclusive B —> X s £ + £~ decay presents a complemen¬ 
tary and also more complex test of the SM since different contributions add to the decay rate, 
see fig. [l|. Quite recently, BELLE announced the first measurement of the inclusive mode [Q . 

Generally, inclusive rare decay modes of the B meson are theoretically clean observables. 
For instance the decay width T(B —> AQy) is well approximated by the partonic decay rate 
r (6 —> sy), which can be analysed in renormalization-group-improved perturbation theory. 
Non-perturbative contributions play only a subdominant role and can be calculated in a model- 
independent way by using the heavy-quark expansion. However, in the decay B —> X s £ + £~ 
there are also on-shell cc resonances. While in the decay B —> X s j (on-shell photon) the inter¬ 
mediate -0 background for example, namely B —«■ ipX s followed by -0 —> A' 7 , is suppressed and 
can be subtracted from the B —> A s y decay rate, the cc resonances in the decay B —> X s l + l~ 
(off-shell photon) show up as large peaks in the dilepton invariant mass spectrum. These res¬ 
onances can be removed by appropriate kinematic cuts in the invariant mass spectrum. In the 
’perturbative windows’, namely in the low-dilepton mass region s = ( p^+ + pe-) 2 /mb 2 < 0.25 
and also in the high-dilepton mass region with 0.65 < s, theoretical predictions for the invari¬ 
ant mass spectrum are dominated by the purely perturbative contributions, and theoretical 
precision comparable with the one reached in the decay B —> AQy is possible. In the decay 
B —> X s £ + £~ kinematic observables such as the invariant dilepton mass spectrum and the 
forward-backward (FB) asymmetry are particularly attractive, especially for the search for 
physics beyond the SM. These observables are usually normalized by the semilcptonic decay 
rate in order to reduce the uncertainties due to bottom quark mass and CKM angles. The 
normalized dilepton invariant mass spectrum and the FB asymmetry are defined as 

R(s) = ^-r (B - X s £ + £~)/T(B - X c ev), (1.1) 



Afb(s) = 


f 1 , „ d 2 T(B -> X s £+£~) . n . 

x / d cos 0i - - — ---sgn(cos^), 

J -i ds a cos tie 


T(B -> X c eV) 

where 0£ is the angle between £ + and B momenta in the dilepton centre-of-mass frame. 


( 1 . 2 ) 
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The B factories will soon provide statistics and resolution needed for the measurements 
of B —> X s £ + £~ kinematic distributions. Precise theoretical estimates of the SM expectations 
are therefore needed in order to perform new significant tests of flavour physics. 

2 NNLL QCD Corrections 

Within inclusive B decay modes, short-distance QCD corrections lead to a sizeable modifi¬ 
cation of the pure electroweak contribution, generating large logarithms of the form a™{mb) 
log m (mf,/Mheavy), where m < n (with n = 0 , 1 , 2 ,...). A suitable framework to achieve the 
necessary resummations of these large logs is the construction of an effective low-energy theory 
with five quarks, obtained by integrating out the heavy degrees of freedom, 

AC 10 

ft// = At £ Oft) Gift) ■ (2.3) 

V 2 j=i 

Compared with the decay B —> X s j, the effective Hamiltonian (|2.3|) contains the two addi¬ 
tional operators of order 0{a em ), O g and O 10 : 


Oi = 

{s^^T a P L c) ( cYT a P L b ) 

o 2 = 

(sq^Pic) (cry^Pib) 

o 7 = 

e/167r 2 m b (n) {sa^Pnb) F^ u 

Os = 

g s /l6ir 2 m b {n) {sa^T a P R b) G\ 

o 9 = 

e 2 /167r 2 (s7 M P L i>) (T/T) 

010 = 

e 2 /167r 2 (sy^Pife) {(-Yld) 


The four-quark operators O 3 6 are not given explicitly because of their numerically small 
Wilson coefficients. 

Within this framework, QCD corrections are twofold: corrections related to the Wilson 
coefficients, and those related to the matrix elements of the various operators, both evaluated 
at the low-energy scale g ~ m b . As the heavy fields are integrated out, the top-quark-, W-, and 
Z-mass dependence is contained in the initial conditions of the Wilson coefficients, determined 
by a matching procedure between full and effective theory at the high scale. By means of RG 
equations, the C)(/i, M^ escvy ) are then evolved at the low scale. Finally, the QCD corrections 
to the matrix elements of the operators are evaluated at the low scale. 

Because the first large logarithm of the form logarises already without gluons 
due to the mixing of the four-quark operator 0 2 into O g at one loop, the leading logarithms 
sum (LL) and the next-to-leading logarithms sum (NLL) are given by 

LL [a em \og{m b /M)] a™{m b ) log n {m b /M) 

NLL [« em log {m b /M)\ < + 1 (m 6 ) log n {m b /M ). 

The complete NLL contributions to the decay amplitude can be found in [|7|, || . Since the 
LL contribution to the rate turns out to be numerically rather small, NLL terms represent an 
0(1) correction to this observable. On the other hand, since a non-vanishing FB asymmetry 
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Figure 2: Typical diagrams (finite parts) contributing to the matrix element of the operator 
CV '2 at the NNLL level 


is generated by the interference of vector (~ 07 , 9 ) and axial-vector (~ 0 1O ) leptonic currents, 
the LL amplitude leads to a vanishing result and NLL terms represent the lowest non-trivial 
contribution to this observable. 

For these reasons, a computation of NNLL terms in B —> X s £Xt~ is needed if one aims at 
the same numerical accuracy as achieved by the NLL analysis of B —> X s ^. Large parts of 
the latter can be taken over and used in the NNLL calculation of B —> X s £ + £~. However, this 
is not the full story. 

The full computation of initial conditions of the renormalization group equation to NNLL 
precision has been presented in Ref. [|l(|] some time ago - including a confirmation of the 


b —> sy NLL matching results of 03 ■ The inclusion of this NNLL contribution removes the 
large matching scale uncertainty (around 16%) of the NLL calculation of the b —> s£ + £~ decay 
rate. 

Most of the NNLL contributions to the anomalous-dimension matrix can be derived from 
the NLL analysis of b —> sy [9|. The missing entries are estimated to have a small numerical 
influence on the dilcpton mass spectrum [|1|] and do not contribute to the FB asymmetry. 

There are two further important ingredients of the NNLL program which were recently 
calculated, namely the two-loop matrix elements of the four-quark operators 0 1,2 and the 
NNLL bremsstrahlung contributions which will be discussed in the following sections. 

In principle, a complete NNLL calculation of the B —> X s £ + £~ rate would require also 
the calculation of two-loop matrix element of the operator Og. However, its impact to the 
dilepton mass spectrum is also estimated to be very small. Similarly to the missing piece 
of the anomalous-dimension matrix, also this (scale-independent) contribution does not enter 
the FB asymmetry at NNLL accuracy. 


3 Two-loop Matrix Elements of O 


1,2 


Within the B —> calculation at NLL, the two-loop matrix elements of the four-quark 

operator 0 2 for an on-shell photon were calculated in [12|] using Mellin-Barnes techniques. 
This calculation was extended in J5| to the case of an off-shell photon (see fig. ]2|) with the help 
of a double Mellin-Barnes representation which corresponds to a NNLL contribution relevant 
to the decay B —> X s £ + £~. This leads to a double expansion in the dilepton mass s and the 
mass ratio m 2 c /m|. Thus, the validity of these analytical results given in || is restricted to 
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Figure 3: Two gauge-invariant subset, without (up) and with (down) a nontrivial cc threshold. 


small dilepton masses s < 0.25. 

An independent check of these results has been performed by us (T^j. Moreover, our NNLL 
calculation is also valid for high dilepton masses for which the experimental methods have 
much higher efficiency compared to the one at low clilepton masses J4J. 

In our approach (T^j, the following calculational method was used: First all diagrams were 
converted into sums of sun-set type integrals and their mass derivatives, 


d n pd n q 


p 1 


,/n 




■ ■q 


A l j 


{{jp ± k) 2 ± m\)(q 2 ± mDir 2 + m 2 ) 

where r = p + q. The effective masses mf 2 3 and the effective momentum k are polynomial 
functions of physical masses, external kinematics, and Feynman parameters associated with 
the diagrams. Integrations over the Feynman parameters are understood. After the internal 
momenta p and q are integrated over, one can show that the results can be spanned by a 
finite set of ten scalar kernels, which are functions of m\ 2 3 and k 2 in one-dimensional integral 
representations, multiplied by tensors made of the metric and k. The Feynman parameters are 
extended into the complex plane to effect rapid numerical convergence. In principle, possible 
IR singularities have to be isolated and substracted. However, an important simplification 
within this specific calculation is given by the fact that all relevant two-loop diagrams are IR 
finite. 

It is clear that a nontrivial cc threshold behaviour at the NNLL level cannot be reproduced 
by the expansion method used in ||. This easily explains our results of the comparison of 
the two independent calculations, |J and [Hj, within the low-s-region. The Mellin-Barnes 


expansion of the gauge-invariant subset shown in fig. [| are in excellent agreement with our 
numerical results. We found that the expanded results given in || are even valid beyond the 
claimed validity range s < 0.25. In contrast, this is not true for the second gauge-invariant 
subset given in fig. [| because of the nontrivial cc threshold in that case. In [Q it was already 
shown that in the low-dilepton mass region these NNLL contributions reduce the perturbative 
uncertainty (due the low-scale dependence) from ±13% down to ±6.5% and also the central 
value is changed significantly, ~ 14%. 

There is no additional problem due to the charm mass renormalisation scheme ambiguity 
within the decay B —> X s i + i~ because the charm dependence starts already at one-loop in 
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Figure 4: Virtual (up) and real (down) QCD corrections. 


contrast to the case of the decay B —> X s j. The charm dependence itself leads to a ~ 7% uncer¬ 
tainty. These small uncertainties in the inclusive mode should be compared with the ones of the 
correspondung exclusive mode B —> /i*/i + /x~ given in JF|; A BR = (1^, ± 6,^4 ,^ 4 , ±2)%. 


The first dominating error represents the hadronic uncertainty due to the formfaetors. 

A phenomenolgical NNLL analysis including the high dilepton mass region will be presented 

0 - 


111 


4 Bremsstrahlung Contributions 

The NNLL bremsstrahlung contributions were also recently calculated for the dilepton mass 
spectrum (symmetric part) in | 6 ], 14] and for the FB asymmetry in [14. 16], see fig. [|. 

In |L| we have separated the bremsstrahlung corrections into universal and nonuniversal 
pieces. We defined the terms cr 7 i 9 (s) which take into account universal 0(a s ) bremsstrahlung 
corrections as 0(a s ) corrections of the effective Wilson coefficients by demanding that the re¬ 
maining (finite) non-universal bremsstrahlung corrections of the rate (encoded in the functions 
T77 and 199) vanish in the limit s —> 1. Then the universal contributions take into account 
the truly soft component of the radiation, which diverges at the s —> 1 boundary of the phase 
space. This is because in the s —> 1 limit only the soft component of the radiation survives and, 
according to Low’s theorem, the latter gives rise to a correction proportional to the tree-level 
matrix element. In |l4j we then found that all nonuniversal bremsstrahlung corrections r,;, to 


the rate and to the FB asymmetry, are rather small all over the phase space, and particularly 
for large values of s (|vy(s)| < 0.5 for s > 0.3), in comparison with the dominating universal 
corrections. 

In the case of the forward-backward asymmetry there are ambiguities arising for d 7 ^ 4 in 
the definition of 75 but in the case of the decay rate, the problematic 75 contribution vanishes 
because of the pi p 2 permutation symmetry of the leptonic phase space. To circumvent 
this problem, we employed the following hybrid regularization scheme |14|: the Dirac algebra 
of IR-divergent pieces is strictly treated in four dimensions (dimensional reduction), while the 
virtual UV-divergent pieces, which do not involve any 75 ambiguity, are still computed in naive 
dimensional regularization. For example, the hybrid wave function renormalization constant 
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Figure 5: Comparison between NNLL and NLL results for Afb(s) in the low s region. The 
three thick (red) lines are the NNLL predictions for fi = 5 GeV (full), and /j = 2.5 and 10 
GeV (dashed); the dotted (blue) curves are the corresponding NLL results. All curves for 
m c /mb = 0.29. 


for a massless quark is given by Z™ =0 = 1 — ^ — l). At this level of the pertur¬ 

bative expansion, this hybrid regularization scheme is gauge invariant. Using this scheme we 
were able to explicitly verify the cancellation of IR divergences in our bremsstrahlung calcu¬ 
lation |b|. 


However, anticipating the cancellation of IR-divergences, one also can circumvent 


the 75 ambiguity by calculating only finite bremsstrahlung contributions [|T 6 ] . 

Let us finally discuss the phenomenological impact of these bremsstrahlung calculations - 
focusing on the position of the zero of the FB asymmetry. This quantity, defined by Afb(so) = 
0, is particularly interesting to determine relative sign and magnitude of the Wilson coefficients 
C 7 and Cg and it is therefore extremely sensitive to possible new physics effects. 


The NLL result 


„NLL 


= 0.14 ±0.02 where the error is determined by the scale dependence 


(2.5 GeV < /j, < 10 GeV) is now modified by the discussed NNLL contributions to (see fig. ||) 


jwll = o ^ 62 ± o ^ og 


(4.5) 


In this case the variation of the result induced by the scale dependence is accidentally very 
small (about ±1% for 2.5 GeV < fi < 10 GeV) and cannot be regarded as a good estimate of 
missing higher-order effects. Taking into account the separate scale variation of both Wilson 
coefficients Cg ew and Cj ew , and the charm-mass dependence, we estimate a conservative overall 
error on so of about 5% |14|. In this s-region the nonperturbative l/ml and 1/m;? corrections 
to A F b(s) are very small and also under control mm PI P|| . Summing up, the zero of the 
FB asymmetry in the inclusive mode turns out to be one of the most sensitive tests for new 
physics beyond the SM. 
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